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Abstract 

We describe various aspects of statistical mechanics defined in the complex temper- 
ature or coupling-constant plane. Using exactly solvable models, we analyse such 
aspects as renormalization group fiows in the complex plane, the distribution of 
partition function zeros, and the question of new coupling-constant symmetries of 
complex-plane spin models. The double-scaling form of matrix models is shown to 
be exactly equivalent to finite-size scaling of 2-dimensional spin systems. This is 
used to show that the string susceptibility exponents derived from matrix models 
can be obtained numerically with very high accuracy from the scaling of finite- 
partition function zeros in the complex plane. 
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1 Introduction 

In contrast with wliat miglit be expected at first siglit, tlie generalization of spin 
models and gauge theories to the complex activity or temperature p| plane is a subject 
rich in physical content. In general, information concerning the phase structure, and, 
in particular, the possible phase transitions for physical values of the parameters can 
be extracted from the behaviour of the finite volume partition function in the complex 
plane. The Yang-Lee edge singularity ^ is one prime example of this, the predictions 
of the distribution of partition function zeros in the complex plane another |^ . That these 
beautiful results can also be very much of practical use in the study of phase transitions, 
has been amply demonstrated 

The purpose of this paper is to consider some aspects of spin, matrix and gauge 
theories in the complex plane that may have gone unnoticed until now. In so far as it 
is possible, we shall compare general predictions for the complex plane behaviour with 
exact analytical solutions. By doing so, we unravel various subtleties involved in defining 
these models for complex parameters, and we also gain additional insight into the nature 
of statistical mechanics off the real axis. 

The simplest example of an exactly solvable model in statistical mechanics is the 
one-dimensional Ising model. Although this model cannot have a real phase transition 
at non-zero temperature, it does in fact have a continuous phase transition at T = 0. 
Because of the special nature of a phase transition occurring at the origin, not all of the 
formalism of critical phenomena can directly be applied to this system. In particular, 
only certain ratios of critical exponents are uniquely defined. Nevertheless, this is often 
sufficient to test the more general predictions of complex plane partition functions. In 
section 0, we use this model to illustrate some of these more general results concerning 
the location of partition function zeros [Q, and use it to comment on some interesting 
aspects of renormalization group fiows in the complex plane. This is possible, since one can 
construct exact renormalization group transformations in the space of just three operators 
(the energy, the spin and the constant operator) by means of spin decimation. 

In section ^ we turn to the question of new symmetries of discrete spin models 
at complex temperature, an issue recently brought up in the case of the Ising model by 
Marchesini and Shrock [0. We show here that the non-trivial part of these new symmetries 
requires a careful definition of what is meant by the thermodynamic limit in the complex 
plane. 

Section ^ contains a discussion of the double scaling limit in matrix models and a 
discussion of 2D U (N) lattice gauge theory in the complex coupling constant plane. The 
parameter A^, which eventually is sent to infinity in order to obtain a phase transition 0, 
is shown to play a role quite analogous to the volume in ordinary systems with a finite 
number of degrees of freedom per site (or plaquette). In this manner, we demonstrate that 
the double scaling of matrix models (and hence low- dimensional string theory) is nothing 
but ordinary finite-size scaling of models in statistical mechanics. This allows us to extract 
numerically the "string susceptibility" exponent of matrix models through the scaling of 
their corresponding partition function zeros. Finally, section |^ contains our conclusions. 
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2 The ID Ising Model at Complex Activity and Complex Temperature 

Define the partition function for a finite number of sites as 



Zn = 



with periodic boundary conditions, i.e., (Jn+i = Introducing the symmetric transfer 
matrix V{a, cr'), it is trivial to compute the partition function from the trace: 

Zn = = A^' + A^, (2) 

where A+ and A_ are the two eigenvalues. Noting that V{a, a') can be represented by a 
2x2 matrix, 



V{a,a') = e 



(3) 
one has 

A± = e^ cosh(/i) ± v'e2^ sinh^ {h) + q-^^ . (4) 

This is all well known. On the real axis A^ > A_, and the thermodynamic limit of the 
partition function, reached as N oo, involves only A+. For example, the free energy per 
spin is found to be 

/ = lim A^^MuZtv = lnA+. (5) 
Similarly, the magnetization per spin is explicitly given by 

m{h,K) = lim iV-iZ^i^ [Zjv] 

N—>-oo on 

sinh(/i) 
[e2^ sinh2(/i) +e-2^]V2- 

It is also well-known that this system cannot have a phase transition for non-zero 
temperatures, and any real value of the magnetic field. We see this explicitly from the 
solution (^). Nevertheless, at h = the system does have a continuous phase transition 
at T = {K = oo), which to a very large extent can be phrased within the standard 
formalism of critical phenomena Since the transition occurs at an asymmetric point 
at which the usual scaled temperature variable t = (T — Tc)/Tc cannot be defined, one 
must be more careful and rewrite all scaling relations in terms of the correlation length. 
This in turn is computed from the 2-point function |^, 

g,, ^ {a,a,)-{a,){a,) = sm\2<P)[X^/ X^f-^\ , (7) 

where is determined from the solution < < 7r/2 of cot(20) = e2^sinh(/i). This 
shows that 

e = l/ln[A+/A_], (8) 

which at /i = is well-defined (and non-singular) for all positive temperatures. The usual 
critical exponents are then defined through the behaviour of the singular parts of the 
thermodynamic quantities as one approaches the critical point at T = 0: 

A/(0,T)~r^'-"^/% m(0,T)~r^/% x(0,T)~e^ (9) 
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From the exact solution (j^) it follows that 



7 = z/ = 2-a, (10) 

while from the corresponding exact solution of the 2-point function we see that rj = 1. 
Finally, from the behaviour of the magnetization at non-vanishing magnetic field at the 
critical point, m{h,0) = sgn{h), it follows that 6 = oo. Although this model does not 
display spontaneous magnetization, it follows from the definition (P) that (3 = 0. These 
values of the critical exponents are all consistent with the general scaling relations 

a + 2/? + 7 = 2 

(36 = /3 + 7 
7 = (2-r/)z/ (11) 

for (35 = u [a result which also follows formally from the exact solution given above), and 
with the general hyperscaling relation 

du = 2 — a 

6 = {d + 2-r])/{d-2 + 7]) (12) 

for d = 1. Thus, although not all critical exponents are uniquely defined in this case, 
certain ratios of exponents are completely well defined, and in full accord with general 
scaling relations at a critical point. 

These simple results do not hold in general, however, once we go off the real axis in 
either magnetic field h (or activity e~^^) or temperature (proportional to K~^). The reason 
is, of course, that the two transfer matrix eigenvalues X± become in general complex- 
valued. The simple thermodynamic limit N ^ oo then obviously does not always exist 
(if, e.g., |A±| = 1, the naive thermodynamic limit of the partition function is given by a 
non-convergent phase running around on the unit circle) . But there are clearly also many 
regions in the complex plane where this naive thermodynamic limit is well-defined. The 
location of partition function zeros Zn = is one such example, which we turn to next. 

One of the predictions of ref. is that for a finite system of volume (number of 
sites) N, the partition function zeros close to the critical point should scale as 

hi = iV-2WVn(0), (13) 

where the index n indicates the nth partition function zero. Also, the scaling function 
fniKN^/") evaluated at zero argument has been predicted to behave, for n large, like @] 

/„(0) = -Cn''"^'^'\ (14) 

where C is some constant. In other words, close to the critical point we should have 

/,2 = _cN-^P&/-n^ms+i)^ (15) 

which in the case of the ID Ising model unambiguously implies 

K = ^C'"^. (16) 
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Let us compare this with the known exact solution. In the hmit K >> 1, it follows 
from eq. (^) that 

A+ = e^cosh(/i) +e^sinh(/i) 

A_ = e^cosh(/i) -e^sinh(/i). (17) 
Thus, in this limit we see from eq. (Q) that 

Zjv ~ 2e^^ cosh{Nh). (18) 
The partition function zeros are therefore to be found at 

or, since n should be taken large in order to compare with eq. (0), 

Tt 

K ~ (20) 



N- 



which is precisely the expected behaviour (|16D, with C = ti 



^2 



In the 2D Ising model the partition function zeros at zero magnetic field lie on 
two intersecting circles in the complex temperature plane . One of these circles crosses 
the real temperature axis at the 2D Ising critical temperature Tc at a right angle. For 
the ID Ising model it follows from eq. @) that the zero-field partition function zeros lie 
on just one circle, the unit circle in the complex tanh(if)-plane, crossing the real axis at 
tanh(i^) = 1 (corresponding to the ID Ising critical coupling K = oo) again at a right 
angle. Although not rigorously related to the general prediction for this angle (p of ref. 

, , cos(7ra) — A_M+ 
tan 2-aV = .\ / , 21 

this behaviour is at least consistent with it. Here a is the usual critical exponent, and A± 



are the specific heat amplitudes above and below T^. For ip = 7r/2 eq. (|2l| ) can be brought 
to 

[1 - A_/A+] cos(a7r) = -[1 - A_/A+]. (22) 

Solutions are either a = 1 with A-/A+ undetermined or A^/A^ = 1 with a undetermined. 
The latter agrees with what we previously found and hence the prediction (|21| ) applies to 
the ID Ising model. 



2.1 Complex-plane Singularities of the ID Ising Model 

Let us define a " critical point" anywhere in the complex plane by the condition of 
a diverging correlation length, = oo. In the ID Ising model we have already seen that 
^ = 1/ ln[A-(./A_], where A-t are the two transfer matrix eigenvalues. The condition ^ = oo 
in this case implies A+ = A_. When h = 0, this gives 

cosh(ir) = sinh(ir), (23) 

which for real K can only be satisfied asymptotically at = oo. This is just the usual 
ID Ising transition at T = 0. 
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When h ^ 0, there are more possibihties. In particular, A+ = A_ in the whole 
subspace where 

e^^ smh\h) + e-^^ = 0. (24) 

Consider first a simple special case, where K is real, but h is allowed to become complex. 
Then, parametrizing h = hi + ih2, and using some trigonometric identities, we find that 
vanishing of the imaginary part of eq. ( p4D implies 

2e2^sinh(/ii)cos(/i2)cosh(/ii)sin(/i2) = 0, (25) 

i.e., hi = 0, or /12 = rnc /2 or both. Demanding also the vanishing of the real part of 
eq. ( plD thus splits up in two parts: (a) hi = and (b) /12 = rni /2. Case (a) implies 

sin2(/i2) = exp[-4i^'], (26) 

which can only be satisfied for i^' > 0. The solution to this equation in the complex 
e'^-plane is shown graphically in fig. |I|a. 
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Figure 1: Critical points for the ID Ising model with complex couplings. 



Case (b) leads to 

e^^ Un\i\hi) cos' (^] - cosh2(/ii 



1 - cos^ ( — 



} + e-2^ = 0, 



(27) 



which in turn has two possible solutions. In the first case, cos^(n7r/2) = +1, for which the 

0. This cannot be satisfied for real hi 
(2j + 1), j = 0, 1, . . ., which leads to 



above equation reduces to e sinh (hi) + e 
and K. The other case is cos^(n'7r/2) 



-2K 



0, i.e. n 



cosh^(/ii 



exp[-4K]. 



f28) 
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Although this equation cannot be satisfied for > 0, there are non-trivial solutions for 
K < 0, which we can think of as the antiferromagnetic situation. In the complex e^-plane 
this corresponds to the solutions shown in fig. |I]b. 

At these points in the complex plane we have not only ^ = 00, but also magneti- 
zation m = 00. The finite volume partition function in these cases equals 

Zn = 2[e^'cosh(/i)]^, (29) 

so that, on account of cosh(/ii + 2/12) = cosh(/ii) cos(/i2) + ismh.{hi) sin(/i2), the case (a) 
corresponds to Zn real. In contrast, in the antiferromagnetic case of (b), A+ = A_ are 
purely imaginary. The simplest thermodynamic limit N ^ 00 does not exist in this case. 

On viewing the above results, one question that immediately arises is: how can a 
local quantity such as the magnetization per spin diverge[^ at these complex-plane critical 
points? 

To answer this question in detail, let us try to regularize everything by considering 
first a finite volume of spins. We still have 

m{h,K) ^ Z^iA/^l (30) 



dh I A^ 

but now both transfer matrix eigenvalues must be kept in the expression for the partition 
function, i.e. = X^ + X^. Carrying out the differentiation, a convenient way to express 
the finite-volume magnetization is as 



e^ sinh(/i) X^-X^ 
[e2^ sinh2(/i) +e-2^]V2 X^TI^ 



Ordinarily, on the real axes, A_|_ > A_, in which case the above expression reduces 
in the thermodynamic limit A^ — > 00 to the magnetization m{h,K) of eq. (^. But the 
present case corresponds to A+ = A_, so we must be more careful. If we take the limit 
A_ A_|_ (where e^^ sinh^(/i) + e~^-^ — > 0) for any finite N, then 

_ ^ [e^ cosh{h)f~' x2Nx (e^^ smh\h) + (32) 

while of course 

A^' + A^ ^ 2[e^cosh(/i)]^ (33) 
to lowest order. Thus at any finite A^ we have 

m{h, K) = ^f^'''^^^^ = N tanh(/i) (34) 
e-^ cosh(/;,) 

at these complex-plane critical points. This behaviour, linear growth with the size of the 
system, is actually precisely as expected from finite-size scaling theory when we treat 

1) The absolute value of the magnetization (P) is clearly bound by unity for all real 
values of K and h, in agreement with the notion that we are making a statistical 
average over a number fiuctuating between +1 and —1. The magnetization escapes 
this bound in the complex plane by being the average over complex- valued Boltzmann 
factor "probabilities" . 
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these complex-plane singularities as true critical points. To see this, consider again the 
infinite- volume magnetization at, for simplicity, purely imaginary external magnetic field 
h = ih2, 



m(/i2, K) 



ie sin(/i2 



(35) 



(36) 





[A+1 




i = ln 


= In 




A_ 





(37) 



-sin2(/i2)e2^]V2 " 

Expanding /i2 = + ^^^2 around a critical magnetic field /i2 given by sin^(/i2 
exp[— 47^"], we note that 

m{Ah2,K) ~ A/i^ 

as A/i2 — > 0, with p = —1/2. 

Similarly, for the correlation length we have 

" e^cos(/i^) +cv/M; ] 
_e^cos{h*2) -cVAh^l ^ 

with u = 1/2. Finite-size scaling theory then predicts 

m{Ah2,K) ~ N-P/" ~ N, (38) 
in agreement with the explicit solution (^4]). 

2.2 Renormalization Group Flows in the Complex Plane 

Not surprisingly, the ID Ising model can be solved by an exact real-space renor- 
malization group transformation, spin decimation. This hence gives us the possibility of 
explicitly observing some exact results concerning RG fiows in the complex plane, results 
that cannot easily be obtained from more general principles. For a rescaling factor b = 2, 
the exact RG transformation for the ID Ising model takes the following form M: 



exp[2h'] 



cosh(2is: + h) cosh(2ir - h) 

cosh^(/i) 
2h cosh(2/r + h) 
^ cosh(2/s: - h) ' 



(39) 



where primed variables are the renormalized couplings. The RG fiow in the K — h plane 
is not trivial [Q, but when h = Q it becomes very simple: we get exp[2/i'] = 1, i.e. h' = 0. 
This just expresses consistency; when we start with no magnetic field, we do not generate 
it in the process of renormalization. The fiow in the temperature direction is slightly less 
trivial: 



cosh2(2ir). 



(40) 



This equation has two fixed points, K* = and K* = oo. The fixed point at K* = cxd is 
unstable, and it is the usual ID Ising critical point; here ^ = oo. 

Now, from general principles the infinite correlation length condition ^ = oo, which 
is fulfilled on the surface 

smh\h) = -e"^^, (41) 

must be associated with a "critical hypersurface" in the complex plane. Once on this 
surface, the RG fiow must remain on it. It is easy to verify that 

e~^^ = — sinh^(/i') if and only if 
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-AK 



sinh^(/i). 



(42) 



Note that this complex-plane hypersurface intersects the plane with real couplings K and 
h at the usual critical point {K,h) = {oo,0), so we should expect that this surface in 
almost all respects can be treated as a standard critical hypersurface. (There should then 
also be complex-plane equivalents of irrelevant operators, pushing the RG flow towards a 
fixed point on this surface). 

It follows from eqs. ( ^9|) that the critical surface given by (|4T|) contains, besides the 
unstable fixed point {K,h) = (oo,0), two stable ones at {K,h) = (0,±i7r/2) (everything 
is symmetric under h ~h). The critical line in the interesting coupling plane (e~^, ih2) 
and the RG flow on it are shown in fig. ^. The figure also shows the results of iterating 
the RG equations starting at several points a little off of the critical line. 
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Figure 2: Renormalization group flow for the ID Ising model in the plane and h = ih2. 
The critical line is shown with the direction of flow indicated. Some results from iterating 
the RG away from the critical line are also given with different plotting symbols and 
connected with dotted lines to guide the eye. The flow is always away from the origin. 

What about the zeros of Zn7 Let us start with a lattice of couplings {K, h = 0) 
and a number of sites equal to 2N, and then apply the RG transformation of rescaling 
factor 2. The zeros of Z2N sue given by 



-K 



jTr{2n+l)/2N 



+ e- 



-K 



(43) 



and squaring this, we get successively 



^i7r(2n+l)/Af 



Ve^ + e-^ ) 
cosh(2ir) - 1 

cosh(2ir) + 1 
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cosh^/^(2fs:) - cosh"^/2(2A') 
cosh^/2(2A:) + cosh-^/2(2A:) 



e 



(44) 



The zeros of the partition function thus behave as expected under this exact RG trans- 
formation: the zeros of Z2N move to those of ^at. 

The RG flow in the complex plane is therefore far from smooth, and except for 
the neighbourhood of the critical hypersurface discussed above, it does not much resem- 
ble the usual RG flows for real couplings. Consider plotting the flow in the complex 
tanh(ii')-plane. On the unit circle of this plane the partition function zeros jump around 
corresponding to going from exp[i7r(2n + 1)/2A^] to exp[2'7r(2n + 1)/A^] as -ft" goes into K' . 
But in the ordinary coupling constant plane of this means going from, for example in 
the last case N = 1: K = ±i7r/4 to K' = —00 in one RG step. 



2.3 The Generalized Heisenberg Model 
The partition function 



Z 



N 



E 



3 ^.7 = 1 • 



Ell ^-"^i 



(45) 



defines a natural generalization of the ID Ising model, - the generalized Heisenberg model 
in ID. Here (fi denotes an n-dimensional vector of unit length, and h is a vector of fixed 
orientation (which can be chosen at will) describing an external magnetic field. For n = 1 
this reduces to the ordinary ID Ising model, while the cases n = 2 and n = 3 correspond 
to the planar and original classical Heisenberg model, respectively. The model is solvable 
for all n in one dimension when the external field vanishes UTTI. We shall here use this 



exact solution to reinvestigate complex-temperature partition functions. 



The eigenvalues of the transfer matrix of eq. (^5]) are known in closed analytical 



form [11, 12 



(46) 



where i = 0, 1, . . ., and is the modified Bessel function of order z/. The model 

can be studied straightforwardly for integer values of n, in which case it yields only 
small modifications in comparison with the Ising case. (The phase transition temperature 
remains at Tc = 0, the critical indices being equal to the Ising indices for all finite values of 
n.) What is more interesting for our purpose here is that the model also admits a unique 



analytic continuation in n for < n < 1 |T2[. In this case the phase transition takes on 
a more familiar form, occurring now at a finite value of T^. The correlation function is 
determined by the two leading transfer matrix eigenvalues, and reads explicitly 

(47) 



This corresponds to a correlation length 



In 



(4^ 
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which diverges at the phase transition point Kc determined by 

I^„{K,) = (49) 

As mentioned above, this phase transition occurs at a finite value of when < n < 1. 
The condition of a diverging correlation length, 

= (50) 

continued into the complex plane thus reduces as n 1 to the circle \coth{K)\ = 1 of the 
ID Ising model. 

Having solvable models with a finite critical point K^, we can investigate the finite- 
size scaling of the zeros of the partition function in the complex plane [Q. The zero 
closest to the real axis should approach the critical point Kc at a rate determined by the 
correlation length exponent i> 

Im(Xo) oc N-^/" . (51) 

As an example we consider n = 0.5 with the critical point at Kc = 0.6755198. The scaling 
(|5l]) is shown in fig. ^. From two successive sizes, we can get an estimate z^(iV) of u. These 
estimates are given in table |l|. 




10 20 50 100 

N 

Figure 3: Finite size scaling of the complex partition function zero closest to the real axis 
for the n = 0.5 ID Heisenberg model. 



We can see that the estimates for u are already quite close to the known exact value 
i> = 1 ||12|. Assuming a finite-size correction of the form A^~^, which fits the data best. 



we can get infinite volume estimates from two successive estimates jy{N). These estimates 
for u, also given in table are remarkably good. 

3 New Symmetries of Spin Models 

In a very interesting paper, Marchesini and Shrock 0, have recently derived what 
is apparently new symmetries of the zero-field Ising model in any number of dimensions. 
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Table 1: Estimates of v from the finite scaling of the imaginary part of the partition 
function zeros computing from successive N (second column) and extrapolated to N = oo 
from two successive estimates (third column). 
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The symmetries are related to discrete transformations of the Ising coupling constant 
(or temperature) in the complex temperature plane. Evidence was found that these new 
symmetries are not without influence on results derived for real coupling constants, since 
they also signal new singularities in the complex plane, singularities that may affect e.g. 
Pade approximants of truncated series expansions on the real coupling constant axis. The 
symmetries, although clearly not of universal nature, are generalizable to other discrete 
spin models. Again, let us consider these general results in the light of the simple solvable 
ID Ising model. We shall later make some related remarks on the 2D Ising case. 

The ID Ising model is of coordination number q = 2, and the two main theorems 
of ref. are therefore applicable in this case. These state that the zero-field {h = 0) free 
energy f{K) should be invariant under the following transformations: 

f{K) = f{K + in7i), n integer (52) 

f{K) = f (^K + i(n + ^^n^ , n integer. (53) 

But if we take the explicit expression for the free energy per spin of the ID Ising model, 

f{K) = -ln[2cosh(ir)], (54) 

this does not appear to be satisfied. In fact, the only invariance of this kind in eq. (Q) is 

K —>■ k + 27rm , n integer, (55) 

which is trivially satisfied by the Ising model in any number of dimensions, on account 
of the integer-valued nature of the spin-spin interaction. How can we understand this 
discrepancy? 
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The explanation lies, of course, in the subtlety involved in defining the thermody- 
namic limit of the partition function off the real axis. In defining the free energy of the 
infinite volume Ising model as in eq. ( |5^ ) above, we have already chosen one particular 
limit, the one valid on the real coupling constant axis. Although this expression can be 
analytically continued into the complex plane, the resulting free energy may not be deriv- 
able from the infinite volume limit of the general partition function in the complex plane. 
As we have already seen in several examples before, these two procedures do not always 
commute. 

To understand the new symmetry, let us therefore revert to the finite volume par- 
tition function Zn = + With zero magnetic field, 

A+ = 2cosh(ir) , A_ = 2sinh(X), (56) 

and consider now the transformation (|53D. Under this, 

cosh(ir) ^ (-l)"isinh(ii') 

sinh(ir) ^ (-l)^cosh(ir), (57) 

or in other words, 

A+ ^ (-1)^A_ , A_ ^ (-1)^A+. (58) 

Consider now taking the thermodynamic limit in steps of 4iV instead of just in 
steps of N , as it is usually defined. Then, with = 4A^, the transformation ( ^3l) acts on 
the partition function as 

[(-1)^A_]^^ + [(-1)"«A+]^^ 

Ai^ + Af = Z^, (59) 

which shows explicitly that this indeed is a symmetry; the transformation has not done 
more than what effectively amounts to swapping eigenvalues. 

But here another subtlety in connection with this new symmetry becomes apparent: 
it is not sufficient to remain at the level of finite volume partition functions, perform the 
transformation, and then take the thermodynamic limit. Instead, the symmetry hinges in 
an essential way on the precise definition of boundary conditions and the very manner 
in which this thermodynamic hmit is defined. On the real axis these problems usually 
do not arise; one obtains a unique thermodynamic limit independently of the manner in 
which the thermodynamic limit is reached.]^ In the complex plane the final thermody- 
namic quantities depend on the way in which the infinite volume limit is taken. From a 
physical point of view this is not acceptable: we must obtain the same bulk thermody- 
namic quantities no matter how we take the infinite volume limit. But again, this problem 
fortunately only surfaces once we enter unphysical parameter values! The need to take 
the thermodynamic limit in units of 4N instead of A^ was noticed already in ref. 0, and 
in fact (as expected from the example shown above), plays a highly non-trivial role in the 

2) Special care is required in the case of spontaneous symmetry breaking, even on the 
real temperature axis. This is perhaps the closest we can get to an ordinary statistical 
mechanics analogue of the complex-T ambiguity discussed above. 
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proof. The problem involved in taking the simplest infinite volume limit ^ oo off the 
real axis becomes apparent if we consider the coupling constant values K + i{n + l/2)7r 
discussed above. If we just let N — »• oo, the partition function becomes a non-convergent 
oscillatory number proportional to e^^'^l'^. 

The same remarks hold for other quantities, but not always with the same re- 
striction on only considering a particular thermodynamic limit. For example, although 
the correlation length ^ = l/ln[A_|_/A_] does change under the transformation (|53D , the 
real part of its analytic continuation into the complex plane remains invariant. Thus, 
= not only at the usual critical point (s), but also at all points in the complex 
plane connected through the transformation (^3]). It was noted in ref. that for the 
2D Ising model the condition Re(^~^) = is fulfilled on the same points in the complex 
plane where the partition function vanishes. This is actually no accident, occurring be- 
cause ^ oc l/ln[Ai/Ao], where Ai^o are the two leading eigenvalues of the 2D Ising transfer 
matrix [|13|. Since the 2D Ising partition function also acquires its zeros when |Ao| = |Ai|, 
this is the rationale for the equality of these two sets. This particular mechanism, degen- 
eracy of two (or more) leading transfer matrix eigenvalues, behind both partition function 
zeros and the condition Re(^~^) = 0, is in fact more general [0. It occurs in the ID Ising 
model as well. Here, without a magnetic field, the condition Re(.^~^) = imphes 

|coth(ir)| = 1, (60) 

while the zeros of are found at 

A+/A_ = e'-(2n+i)/A^^ (61) 

i.e., for the N ^ oo limit again given by the unit circle |coth(i^)| = 1. In fact, it is 
precisely because it is also a partition function zero that the condition Re(^^^) = in the 
complex plane can be associated with an unambiguous thermodynamic limit. 

Similarly, the transformation (|52|) also only defines a real symmetry of the ID Ising 
model if a very particular definition of the thermodynamic limit is used. The case n = 
even is trivial, so we need only consider the case n = odd, i.e. : K —>■ K + {2n + l)i7r. 
Under this transformation, A± — — A-t, and the partition function is symmetric if one 
takes the thermodynamic limit in units of 2N . (Because the coordination number is g = 2, 
there is no need to take the thermodynamic limit in multiples of AN; this also follows 
from the general proof in ref. 0). 

The existence of these kinds of symmetries appears to have been known for some 
time by practitioners of transfer matrix techniques (see, e.g, chapter 7 of ref. 0). 

4 The Double-Scaling Limit: Matrix Models 

The non-analytic behaviour around phase transitions always arises as a result of a 
singular procedure, the thermodynamic limit, in which the number of degrees of freedom 
diverges. For a finite number of degrees of freedom, the partition function is an analytic 
of the coupling (or temperature). It is well known that a thermodynamic limit can be 
achieved not only by letting the volume (or number of sites) go to infinity, but also 
in finite volumes (or even on just points, i.e. zero dimensional field theories) provided 
the number of degrees of freedom diverges by some other means. One example of this 
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is large- matrix models in zero dimensions, which can display behaviour completely 
analogous to ordinary infinite- volume phase transitions. Another classical example is two- 
dimensional QCD in the N oo limit. By gauge invariance, two-dimensional latticized 
SU{N) Yang-Mills theory is a one-plaquette theory. This "trivial" theory with Wilson 
action nevertheless undergoes a 3rd order phase transition in the limit A^ = oo [Q. One 
curious aspect of this 3rd order Gross- Witten transition is that the correlation length, 
defined in a natural physical way as 

e ~ (g2) 

where a is the string tension, remains Gnite at the phase transition. This is difficult 
to reconcile with renormalization group ideas that seem to always require some diverging 
correlation length at non-trivial fixed points. Does this mean that the Gross- Witten phase 
transition cannot be understood from a renormalization group perspective? Or do we just 
have to search a little more for the true underlying correlation length that diverges at the 
transition? 

The purpose of this section is therefore threefold. First, we shall analyze the matrix 
model phase transitions from the point of view of scaling of partition function zeros. 
Second, we shall answer the question: What is the analogue of the "double scaling limit" 
T5| of such matrix models in ordinary systems in statistical mechanics? And third. 



we shall describe the appropriate renormalization group scenario for the Gross- Witten 
transition (which actually does correspond to a diverging correlation length, to be defined 
shortly) . 

We start with the last two questions, which turn out to be intimately related. But 
first a few definitions: By the double scaling limit of matrix models (and low-dimensional 
string theory) we mean the following. The singular part of the string partition function 
({A, Q) (a somewhat confusing terminology in this context, since it will be identified with 
the free energy of certain matrix models) scales as : 

oo „ 

C(A,6;)~(A*-A)2-^«^aG[(A*-A)6;-2/^^] , (63) 

G=0 

where the sum runs over genus G. Since the same factor appears for all genera, this means 
that 

C(A, g) ~ (A* - A)2-^°F((A* - K)g-^'^^) , (64) 

where the scaling function F is simply the sum on the r.h.s. of eq. (|63D. Here A is the 
cosmological constant, and g is the string coupling constant. The exponents 70 and 71 
are known from continuum calculations UTi 



= 12 



25-C + ^(1 -c)(25-c) 
70 + 71 = 2, (65) 



where c is the central charge. We shall mostly consider the simple case c = 0, which 
corresponds to zero-dimensional matrix models. 

In such a matrix model representation, (A* — A) ~ [g* — g) and g ~ 1/A^, where 
g is the matrix model coupling constant, and is the size of the matrix. The double 
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scaling limit entails sending g ^ g*, the critical point of the matrix model at = oo, 
while simultaneously sending oo in such a manner that the scaling variable 



is kept fixed. 



x=ig- g*)N' 



L-1 



2/71 



(66) 



Figure 4: Illustration of a real-space renormalization group transformation. 



What is the analogue of this double scaling limit in ordinary spin systems? It turns 
out to be useful to think of a system at finite size L, i.e., volume V = L"^ in d dimensions. 
A two-dimensional analogue of such a system is shown in fig. ^. Now let us perform a 
simple real-space renormalization group transformation in the form of a spin decimation. 
Wishing to retain the same geometry (although in rescaled form), we choose to integrate 
out all spins along one row and one column, as also indicated in fig. Let us denote 
the spins we integrate out by (p, and the remaining ones by 0'. This spin decimation 
induces a renormalization group fiow; the original Hamiltonian H of operators Oi and 
corresponding coupling constants gi is mapped into a new Hamiltonian H' with in general 
new interactions and in any case new coupling constants g[. That is. 



-H'(<l>',g') 



-H{<t>,<l>',g) 



(67) 



where we have explicitly, as is customary, extracted the constant part G{g) of the new 
Hamiltonian. Defining the free energy by 



f(g)^l-JnZ{g), 



(68) 



this indeed corresponds to the standard real-space renormalization group equation for /, 

f{g) = G{g) + b-'f{g') (69) 

with a length rescaling factor 

(70) 



L-1 

But this time this equation has the meaning of a finite-size renormalization relation. In the 
limit of large volume, L — » oo, the number of spins integrated out at each step compared 
with the number of those retained becomes infinitesimal. This allows us to phrase the 
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real-space renormalization group equation in infinitesimal form, giving in effect a Callan- 
Symanzik-type equation for finite-size scaling. 

The constant part G{g) of the Hamiltonian is assumed to be regular. We write it 

as 

G{g) = \n[Xig)] (71) 
where X{g) must approach unity for large L: 

X{g) = 1 + ^r{g) + o(^). (72) 



Similarly, for large L, the coupling constants g'^ will be close to the previous values gi, so 
that we may define functions j3i{gi) by 

9'i = 9i + \^3^{g^) + o(^]. (73) 



It is then straightforward to see that the recursion relation 

[\{g)]^'ZiL-l,g') = Z{L,g) (74) 



for the partition function implies the following differential equation for the free energy 
f{L,g): 



fiL,g)=r{L,g). (75) 



This is the Callan-Symanzik equation for finite-size scaling. The /3-functions (3i{gi) have 
been used by Roomany and Wyld |jl8| (see also ref. |l^) to get infinite- volume results 
from finite-size scaling. Their main advantage is that they have zeros even in finite vol- 
umes, where no phase transitions can occur. The zeros of these /3-functions determine the 
location of the phase transitions in the infinite-volume systems. 

The Callan-Symanzik equation (^) above gives an exact differential equation for 
the free energy f{L,g) for all values of the couplings gi in the limit of infinite volume, 
L — > oo. It is not restricted to the critical regime close to a critical point given by {g*}- 

For large but finite volumes L"^, we can extract the singular behaviour of the free 
energy near the infinite- volume fixed point {g*}- For simplicity, assume that we have only 
one coupling constant g. Define an ordinary critical point g* by a solution of 

/3ig*) = , /3'ig*) > 0. (76) 

Then the Callan-Symanzik equation (|75D gives the following solution for the singular part 
of the free energy f{L,g) near g*: 

/(L, g) ~ (g* - gy/^'^^'^Fiig* - g)LP'^^*^) , (77) 

where F is an arbitrary scaling function. This of course coincides with the usual finite-size 
scaling formula if we identify the correlation length exponent as follows: 
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an identification that can also be made explicitly by considering the rescaling of the 
correlation length under the spin decimation described above. Note that the regular part 
r{L,g) of the Callan-Symanzik equation ( [75| ) does not influence the singular part of the 
free energy, as expected. 

While this formulation thus allows us to recover the standard finite-size scaling 
formula, it actually also gives us more. It relates the divergence of the correlation length 
directly to the rescaling of couplings through their /3-functions. This will turn out to be 
useful when we return to matrix models. 

Implicit in the above considerations was the assumption of an analytic expansion 
of the /3-function around g*, 

m = Pig*) - f3'{9*){9* -g) + (-1) (79) 

and we also had to assume that l3'{g*) > 0, since otherwise we would not arrive at 
a positive correlation length exponent. A negative correlation length exponent would be 
inconsistent with the assumption of g* being a critical point. More exotic phase transitions 
are, however, also contained in this formalism. Such phase transitions correspond to more 



unusual behaviour of the /3-functions WE 



Note that the finite-size scaling form of the free energy f{L,g) has an uncanny 
resemblance to the Knizhnik-Polyakov-Zamolodchikov scaling form of the string partition 



function (p^) . For this identification to be made, we see from the second relation of eq. (pSD 
that we have to consider the string (matrix model) case as corresponding to ci = 2, and 
identify the size of the matrix with the size L of the system described above. This is 
actually entirely natural, since as far as scaling in size is concerned, a matrix acts like 
a 2-dimensional system. The precise identification is then 

To conclude, the double scaling of matrix models is precisely the hnite-size scaling of 
ordinary 2-dimensional statistical mechanics. 

In matrix model language, the finite-size scaling - as we have seen, the "double 
scaling" - arises in the limit of N —>■ oo, when one simultaneously tunes the matrix 
coupling g to a critical value g* in such a manner that the scaling variable 

x^{g*- g)N'/^^ (81) 

is kept fixed. The derivation of finite-size scaling presented here for spin systems was 
first done for matrix models by Carlson [IT3 through the entirely analogous procedure of 



integrating out one row and one column of an A^ x A^ model, thereby relating it to an 
(A^ — 1) X (A^ — 1) model with renormalized coupling constants. It has very recently been 



rediscovered by Brezin and Zinn- Justin |^0[) who noted the identification ( pOf ) above (see 



also ref.[0]). In fact, in order to facilitate the comparison with the work of Brezin and 



Zinn- Justin [£0[) we have used a notation that closely parallels theirs in the language of 
matrix models. 

This interpretation of the double- scaling limit of matrix models implies immediately 
an almost trivial relation between the string susceptibility 70 (and, equivalently, 71) and 
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the correlation length exponent v: 



z/=^(2-7o)=7i/2. (82) 

One obvious generalization of these relations is to rank-d tensor models (where 
"tensor" is just a loose terminology for any d-index object). Vector models have already 
been analysed from the double-scaling point of view , and found to behave very simi- 
larly to the matrix models. Tensor models of higher rank tensors have also been suggested 



in relation to simplicial gravity in higher dimensions We see from the straightfor- 
ward scaling arguments above that for simple zero-dimensional tensor models, the double 
scaling limit is uniquely given by the rank of the tensor. It is natural to parametrize the 
critical exponents of rank-ci tensor models slightly differently from the matrix model case, 
in order to avoid factors of d in some of the expressions. Let us therefore write the free 
energy for rank-c? tensor models in the form 

/(AT, g) ~ (/ - gf~^^F{{g* - g)N'"^-) . (83) 

This, then, defines the two exponents 70 and 71. Comparing with eq. ([7^), we note that 
with this definition, the relation 

2 - 70 = 71 (84) 

remains valid for all rank-rf tensor models. The correlation length exponent v is then 
related in a very simple manner to the string susceptibility exponents 70 and 71: 

z/=(2-7o)/rf = 7i/^i. (85) 

The correlation length has a very direct interpretation in terms of the rank-rf tensor 
models. It describes the range of effective correlations within invariant subgroups of the 
large- tensor, and it diverges at the critical point. Without this diverging correlation 
length, the description of the double scaling limit in renormalization group terms would 
presumably be void of meaning. 

It is interesting to compare these simple scaling laws with what has already been 
derived by other means. We have earlier seen the direct relation between finite-size scaling 
of matrices and the Knizhnik-Polyakov-Zamolodchikov formula for strings. Another class 
of models which have been studied in the double- scaling limit is that of vector models 



22| , which in our notation corresponds io d = 1. The critical exponents for these vector 
models are precisely given by the relations ( pOD and (p5| ) with d = 1. 

The interpretation of the correlation length in terms of the original "dual" picture 
(random polymers {d = 1), random surfaces {d = 2), etc.) is unfortunately still rather 
obscure. The divergent correlations occur within the rank-d tensor itself, and does not 
seem to have any simple physical meaning. Similarly, if we return to the solvable case of 
2-dimensional lattice-regularized Yang-Mills theory, there is no contradiction between the 
finite physical correlation length (|6^) which is defined in terms of the lattice correlator, 
and the divergent correlation length occurring among subsets of matrix elements inside 
the SU{N) matrix as A^ ^ 00. The latter has no simple interpretation in terms of the 
usual Yang-Mills potentials. 
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Nevertheless, having a well-defined diverging correlation length available makes it 
possible for us to make use of the full machinery of the theory of complex-plane partition 
function zeros discussed in this paper. This in particular implies that also the double- 
scaling limit itself can be entirely understood from the point of view of singularities in 
the complex coupling constant plane. This is the subject of the next subsection. 

4.1 Two-dimensional QCD in the large-A^ limit 

As an application of the above ideas, we shall now demonstrate how the scaling of 
complex coupling constant singularities can be used to numerically determine, e.g., the 
string susceptibility exponents. The single crucial ingredient is the existence of the rela- 
tion (^) between the correlation length exponent and the matrix model susceptibility 
exponent 70. 

For numerical purposes it turns out that the unitary matrix model - 2-dimensional 
latticized Yang-Mills theory - is the easiest matrix model to work with. The relation 
between finite-size scaling and the double-scaling limit as discussed above was derived 
only in the case of a Hermitian matrix ensemble, where the integration over a row and 
a column does not spoil hermiticity. In the unitary ensemble, the mapping onto itself no 
longer holds. For finite values of N, the size of the group elements in the fundamental 
representation, one is not automatically left with a unitary matrix after integrating out 
a row and a column. On the surface, this would seem to invalidate the direct connection 
with finite-size scaling. Indeed, an exact integration of a row and column can only yield 
a unitary matrix model modulo CP^~^ factors |19|]. Fortunately, these corrections are 



subleading in the large- iV limit, and the above formalism therefore does carry through in 
the limit N 00. 

As we are interested in the large-A^ limit, where the distinction between SU{N) 
and U{N) is immaterial, we restrict ourselves to 2D lattice gauge theories of the U{N) 
kind. With the Wilson action one can write the reduced partition function as 

Zj, = /[d[/p]e^"'[^'-^^+^'^^^'l, (86) 



where variables have already been changed to fundamental plaquettes Up, and where a 
trivial volume factor has been taken out. This partition function can be written in closed 
form for any A^ by means of modified Bessel functions |]2^ : 



Zn = detiMij) = det(/i_, (2/(72)) , i^j = 1,2,..., N, (87) 

and it is this expression we shall use to determine the partition function zeros. In order 
to be able to take a smooth large-A^ limit, we rescale variables and introduce A = g'^N. 
The Gross- Witten phase transition occurs at A* = 2 |^. 

The critical exponent 71 from the double-scaling limit of this unitary matrix model 
was computed soon after the discovery of the double-scaling limit in the Hermitian matrix 
ensemble [^. In the notation of this paper, the result is 71 = 3. This corresponds, 
according to eq. (p2D, to a correlation length exponent 

z/ = 3/2 . (88) 

Next, let us consider the model (^) in the complex A-plane. According to the 
general scaling theory of complex-plane zeros of the partition function 0], the partition 
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function zeros closest to the Gross- Witten transition point A* = 2 approach this transition 
point at a rate determined by the correlation length exponent u:^ 

Im(l/Ao) oc N-^/" . (89) 

One can easily show that 

Z^(A) = Z^(A)=Z^(-A). (90) 

Thus each zero of comes in quadruplets, Aq, — Aq, Aq, and — Aq, unless it is purely 

imaginary. We have computed the zeros 1/ A with the smallest imaginary part for = 1 to 
10. Most of these are already listed in but for completeness we give the representative 
in the first quadrant in table 



Table 2: The partition function zeros of eq. closest to the real axis. 



N 


Re(2/Ao) 


Im(2/Ao) 


1 


0.0 


2.404825 


2 


0.639801 


1.490191 


3 


0.810578 


1.130558 


4 


0.883926 


0.930444 


5 


0.922905 


0.800366 


6 


0.946328 


0.707894 


7 


0.961580 


0.638198 


8 


0.972091 


0.583456 


9 


0.979646 


0.539118 


10 


0.985256 


0.502343 



One can now use the finite size, i.e., in our case finite A^, scaling relation (|89|) to ob- 
tain estimates of u. Such estimates, obtained from and — 1 are given in table |^. These 
estimates still have subleading finite A^ corrections. Assuming a simple 1/A^ behaviour we 
can obtain, from two successive i/(A^)'s estimates for the true, infinite A^ values: 

u{oo) = u{N) + {N - l)[iy{N) - u{N - 1)] . (91) 

These estimates for z/(oo) are also listed in table ^. As one can see, the convergence to 
the theoretical value 3/2, eq. (^8|), is very fast. This can also be seen from fig. ^ 

From the exact solution by 0] we can also see that a = —1. Using = 2 as 
appropriate for the matrix model considered, we can see that a and v satisfy the usual 

3) At the upper critical dimension, this behaviour is modified by logarithmic corrections 

4) The exact analytical form of the domain of partition function zeros in the N = oo 
theory can in principle be calculated analytically using the techniques of ref.[^. This 



is an interesting open problem which remains to be solved. The approximate form of 
the partition function zeros in the N = oo theory has been conjectured in ref. p7 
the basis of numerical evidence for finite values of A^. 



on 
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12 5 10 

N 

Figure 5: "Finite N" scaling of the complex partition function zero closest to the real axis 
for the Gross- Witten model. 



(hyper-) scaling relation, eq. ([l^). If this is not fortuitous, and we see no reason why this 
should be the case, it implies that the string susceptibility exponent 71 of the double- 
scaling limit could have been obtained directly from the Gross- Witten solution, using the 
above identifications. 

An interesting question concerns the possible existence of an order parameter de- 
scribing the phase transition. With the correct notion of a magnetic field operator, one 
could establish the missing critical index, and derive the corresponding scaling relations. 
It has been suggested that the Gross- Witten phase transition can be associated with 
the symmetry breaking of finite U{N) subgroups of the full N ^ 00 theory. However, 
the "symmetry breaking" of ref. ||29| occurs on both sides of the phase transition point, 
and it is not immediately obvious how it connects to the conventional notion of an order 
parameter. This aspect of the N = 00 phase transition remains to be better understood. 



5 Conclusions 

The generalization of spin, matrix and gauge models to complex coupling constants 
yields new and useful insight into the dynamics of these systems. We have discussed 
a selection of solvable or near-solvable models where predictions of the complex-plane 
behaviour can be tested directly. The one- dimensional Ising and generalized Heisenberg 
models are particularly useful examples of theories with non-trivial phase transitions 
(albeit, in the Ising model case, at zero temperature) where one can study the behaviour 
in both the complex activity and temperature planes. By exact decimation, one can also 
describe the behavior of renormalization group trajectories in these complex planes. 

We have discussed some of the subtleties involved in classifying new symmetries of 
such spin models in the complex temperature plane. The question of whether these models 
have new symmetries hinges on the precise definition of what is meant by the thermody- 
namic limit for theories with complex Boltzmann factors or transfer-matrix eigenvalues. 
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Table 3: Estimates of from the "finite N" scaling of the imaginary part of the partition 
function zeros computing from successive N (second column) and extrapolated to N = oo 
from two successive estimates (third column). 



Matrix models have recently been under intense study in connection with the 
"double-scaling" limit of low- dimensional string theory. We have clarified various aspects 
of this double-scaling limit by showing that it is equivalent to what in standard statistical 
mechanics terminology is known as finite-size scaling. Here "size" {i.e. volume) is nothing 
but the size of an x matrix. The double scaling of matrix theory is thus equivalent 
to the finite-size scaling of two-dimensional systems, while the corresponding double scal- 
ings of rank-c? tensor models are equivalent to finite size scalings of (i-dimensional systems 
in statistical mechanics. This remarkable connection allows us to phrase the string suscep- 
tibility exponent of matrix models directly in terms of the correlation length exponent u. 
By entering the complex coupling constant plane, the scaling of partition function zeros 
near the critical point is thus given in terms of the string susceptibility exponent. Con- 
versely, we have demonstrated how this scaling of partition function zeros can be used to 
give remarkably accurate values for the string susceptibility exponent. As a by-product, 
we have established that the 3rd order Gross- Witten phase transition of 2-dimensional 
U{N) {N oo) lattice gauge theory with Wilson action is characterized by an underlying 
divergent correlation length, despite the finite value of the naive correlation length 0. 
The Gross- Witten phase transition of 2-dimensional Yang-Mills theory at large A^ can 
thus be understood in entirely conventional renormalization group terms. 



z/(oo) 



2 1.448363 

3 1.468050 

4 1.476775 

5 1.481765 

6 1.485009 

7 1.487288 

8 1.488970 

9 1.490265 

10 1.491294 



1.507424 
1.502950 
1.501725 
1.501227 
1.500966 
1.500745 
1.500620 
1.500561 
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